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I. INTRODUCTION 



In a recent study of the pp — > pp7r° reaction at energies close to threshold, it is found that the angular distributions 
of the outgoing particles are isotropic, in agreement with the assumption that the reaction proceeds as a 33 Po — ► 31 SqSo 
transition |l|J^| . Several model calculations |5|-[| of S-wave pion production, which are based on a single nucleon and 
a pion rescattering mechanism, under estimate the cross section by a factor of 3-5. Inspired by a study of /?-decay 
in nuclei, which indicates that the axial charge of the nuclear system is enhanced by heavy meson exchanges and, by 
using a simple operator form of the NN potential, Lee and Riska [Q have shown that meson exchange currents could 
explain the scale of the cross section. Adding to the single nucleon and rescattering terms a Z-graph describing various 
heavy meson exchanges, Horowitz et al. || performed similar calculations based on explicit one-boson-exchange (OBE) 
model for the NN interaction as well as for the evaluation of meson exchange contributions. 

The contribution from the rescattering term depends on the off-shell behavior of the ttN — ► irN scattering amplitude. 
There are several approaches based on field theoretical models which allow such extension to be made but the results 
are model dependent. In the traditional phenomenological treatment the off mass-shell and on mass-shell 

amplitudes are of the same order of magnitude. Recently, several groups []Qf— [l3f| have shown that the off mass-shell 
7rN rescattering amplitude is enhanced strongly with respect to the on mass-shell amplitude. Hernandez and Oset 
H by applying current algebra and PCAC constraints argue that this enhancement may bring the calculated cross 
section for the pp —> ppir reaction into agreement with experiment. However based on detailed momentum-space 
calculations, Hanhart et al. conclude that the enhancement of the rescattering amplitude due to offshellness falls 
a bit too short to explain the scale of the pp — > ppir cross section. Park et al. jyj, Cohen et al. (lj] and Sato et al. 
Jk| have applied a chiral perturbation theory (^PT), including chiral order and 1 Lagrangian terms. They have 
shown that the off mass-shell 7rN elastic scattering amplitude is enhanced considerably but has an opposite sign with 
respect to the on mass-shell amplitude. Because of this difference in sign, the rescattering term and the Born term 
contributing to the pp —> ppir reaction interfere destructively, making the theoretical cross sections much smaller 
than experimental values and thus suggesting a significant role for heavy meson exchanges in ir° production. 

In the present work we report on a covariant OBE model calculations based on the mechanism depicted in Fig. 1, 
where a virtual boson B (B=7r, a, r/, p, ui, S ...) created on one of the incoming nucleons, is converted into a ir° meson 
on the second via a BN — > 7r°iV conversion process. The half off mass-shell amplitudes for the conversion processes, 
hereafter denoted by T BN ^^o N , are taken as the sum of three terms corresponding to the s, u and t-channels displayed 
in Fig. 2. One important aspect of this mechanism is that at the ir° production threshold, the transferred 4-momentum 
is space-like, q 2 — — 3.3/m~ 2 . This is very much the same kinematic as occurring in ir° electroproduction through 
vector meson exchanges, and we apply a formalism similar to the one applied for electroproduction amplitudes 111] 
calculations. The formalism to be applied is consistent with the OBE picture of the nuclear force and accounts for 
relativistic effects, crossing symmetry, energy dependence and nonlocality of the hadronic interactions. 

We start with the assumption that the pion production in pp — > ppir can be described by the graphs of Fig. 3. The 
diagrams 3a and 3b represent nucleon pole terms where the neutral pion is produced on a nucleon line. The diagram 
3c depicts a pion production occurring on an internal meson line. By applying s-channel unitarity to NN elastic 
scattering in the energy region below two-pion production threshold, it can be shown that mechanisms as such should 
contribute to the production process also |1| . ScMtz et al. reported on a and p meson t pole contributions to 7rN 
scattering. Their model was based on correlated 2tt exchange and constrained by using quasiempirical information 
about the NN — * irir amplitudes. Hanhart et al. used the T-matrix obtained in this model to calculate the 
rescattering term for pp — > ppir . More recently, van Kolck et al. [|l7| considered other short range meson exchange 
mechanisms like p — > lott and <5 — > r/7r° . To keep our model calculations consistent with the OBE picture of Machleidt 
fl9f we limit our discussions to contributions from a, 5 and p poles. In variance with previous calculations jl6 17 
we write these vertices in a more general form and determine the relevant coupling constants by applying the Adler's 
consistency conditions and/or fitting to some observables. We shall demonstrate below that a a- meson t pole term 
which accounts ef fectively for isoscalar-scalar two-pion exchanges contributes significantly to tt° production. With 
this contribution included the model explains the cross-section data for the pp — > ppir reaction at threshold. 

The present article is organized as follows. In Sect. II we present details of the formalism and model parameters. 
The evaluation of the T^ajv—^ojy, T r)p ^, 1K a v and T pp ^^o p conversion amplitudes require knowledge of the airir, Srjir 
and pun vertices. These are deduced in subsections B-D. In Sect. Ill we write the S-wave production amplitude in a 
form suitable for numerical calculations. Model predictions are given in Sect. IV. At energies close to threshold, final 
state interactions (FSI) influences the energy dependence of the calculated cross section [g0[ . These and initial state 
interactions (ISI) corrections are introduced in Sect. V where comparison with data is to be made. We conclude and 
summarize in Sect. VI. 
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II. THE MODEL 



We use the following Lagrangian interaction: 



L = ^Ln^^O^tN + l^-N^^d^N + g aNN NaN + 



n 



9 p nnN (y + ^^ v d v ) TppN + g uN NN-fupN + gs NN NTSN , (2.1) 

with obvious notations. This expression includes terms identical to the ones used by Machleidt et al. jl^] to fit NN 
elastic scattering data in the energy region 0-420 MeV. Here as in Ref. [|l9| pseudovector couplings are assumed for 
the pseudoscalar mesons and the u> tensor coupling is taken to be zero. All of the coupling constants, meson masses 
and cut off parameters are taken from Table A. 2 of Ref. |19| ; their potential C parameter set. 

To calculate the transition amplitude we assume that the reaction is dominated by the mechanism depicted in Fig. 
1 and write the primary production amplitude as, 

M^=J2[ T BN^ON{pi,k; P 2,q)G B (q)S BNN (p 3 , Pl )} + [1^2;3^A] . (2.2) 

B 

Here pi, q and k are 4-momenta of the i-th nucleon, the exchanged boson and the outgoing pion. The sum runs over 
all possible B bosons that may contribute to the process. To be consistent with the OBE picture of the NN force [ p)| |, 
we include exchange contributions from all of the 7r, 77, a, p, lu and S meson s. T he bracket [1 <-> 2; 3 <-» 4] stands for 



a similar sum with the pi, p 3 and P2, Pi momenta interchanged. In Eqn. 2.2 , T BN ^^o N represents the conversion 
amplitude corresponding to the BN — > 7r°iV process. The quantities G B (q) and S B nn(P3,Pi) are the propagator 
and source function of the meson exchanged, respectively. Throughout this work we use covariant expressions for the 
meson propagators and form factors as defined in Ref. EJ. The source functions for scalar, pseudoscalar and vector 
mesons are 

Ssnn( P i,P3) = u(p 3 ) I u(p{) F s (q) , (2.3) 
Spnn(pi,P3) = w(p 3 ) 7 5 I u(pi) F P (q) , (2.4) 

S VNN ( P i,P 3 ) = u{p 3 ) [ 7 "4%?3) + i<r^4 2) (9? 3 )] Iu{p x ) , (2.5) 

where u{p) stands for a nucleon Dirac spinor; / is the appropriate isospin operator and p 3 = p% — q is the final nucleon 
momentum. The functions Fs(q) and Fp{q) are source form factors for scalar and pseudoscalar mesons. For vector 

mesons there are two such quantities Fy^ and Fy^ representing vector and tensor form factors, the analogous of the 
nucleon electromagnetic form factors. In the calculations to be presented below all source form factors are taken in 
the form ||, 

A 2 2 

F B (q 2 ) = 9 BNNfB{q); f B (q) = g ~ ™f . (2.6) 

A% - q A 

Albeit, the propagators and source functions are rather well determined from fitting NN scattering data Thus 
to a large extent, the model success in explaining cross section data for the pp — ► ppn° reaction depends on how well 
the conversion amplitudes Tgjv— >7r° n are calculated. The relative importance of the various exchange contributions 
depend upon the off mass shell behavior of these amplitudes. For example, though very small on mass-shell, the 
amplitude for tt° P — > n°p process is strongly enhanced off mass-shell ||[l0) giving rise to a dominant contribution to 
the production rate. 



A. Nucleon pole s and u channel amplitudes 

In this subsection we write expressions for s and u-channel nucleon pole terms (diagrams 2a and 2b). These are 
common to all of the various meson exchanges. We call p and p' the momenta of incoming and outgoing nucleon, q 
and k those of the boson B and the pion produced, respectively. Let s = (p' + k) 2 be the total energy and 

A "(*) = 5 L T ■ • (2-7) 

x z — M z + te 
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Using the vertices of the Lagrangian, Eqn. 1, and the usual Feynman rules the contributions from diagram 2a with 
B designating a scalar (S), a pseudoscalar (P) and a vector meson (V) are respectively, 



rp(S) _ 

1 s — 



ITU 



M 2 



T 



(s) _ ■ IpNnIttNN 



{p' + kf 

pf + # + M 



X, 



p 

W 



V 



9VN N fnN N f ( \ f n \-f >\ 5u M+ft + M 
-_Mff)A(*)«(p)7 ^ + k) 2_M2 



7 P + 



2M 



u(p) 



where </l — a^q^. After some algebra one obtains, 

rp(s) _ 

1 s — 

.gSNNfirNN 



fs(q)U(k)u(p'h 5 [(» ~ M 2 ) + 2MJf\ u(p)A N (s) 



T { p s) = (2M) 



u{p') 



2 f PN N f-nN N 

mprrin 



M 2 



2M 

rp(s) _ . gVNNf-nNN 

ly — I 



fp(q).U(k) 

M 2 + AM 2 
AM 2 

A N (s)f v (q)U(k)u(p') 7 5 



A N (s)u(p), 



M(s - M 2 ) + 4fc^ - 2M^\ 



1 



(—) 

\2MJ 



(2.8) 
(2.9) 
(2.10) 



(2.11) 



(2.12) 



(2.13) 



The contributions from the u-channel (diagram 2b) are obtained easily by replacing the total energy s in Eqns. 11-13 
with u = (p — k) 2 . Likewise, for a vector meson the u-channel pole term is 



X, 



(«) 



i 9VNN J* NN - A N (s)f v (q)U(k)u(p'h 5 r (l - ^i) [u M 2 + 2M#] u(p) . (2.14) 



B. Evaluation of Tlo 



7T JV-.7T JV 



We now turn to calculate the amplitude T^on^^n- in this case in addition to s and u-channel nucleon pole terms, 
a neutral pion can be formed on an internal er meson line. These terms are displayed graphically in Fig. 4. The 
contribution from graph 4c depends upon the airn meson vertex. There is no information about this coupling from 
traditional analyses of NN elastic scattering |ll| . In order to determine this vertex we apply some physics beyond the 
OBE picture of the NN interactions. Quite generally we write an effective airir vertex in the form, 



V a7TV (k,q) = m a goa + 



k 2 



-gio 



(q-k) 2 



-92a 



(2.15) 



where ffoo-, 5i<r, 52o- are (as yet unknown) constants. This is a slightly generalized Weinberg-type low energy expansion 
plf allowing for three rather than one different constants. Taking the appropriate s and u-channel contribution from 
Eqn. 12 and adding the contribution from a a pole in a t-channel one obtains, 



{2MU 1 



tNN 



\ TO. 

Uq-k) 
(q — k) 2 — m 2 



fAq)U(k)u{p') 



i 



i 



M 2 u-M 2 



1 

M 



u(p) 



V avn u(p')u(p) 



(2.16) 



Now, in order to determine the vertex constants go a ,gi a and g^a we require that the conversion amplitude obeys the 
three Adler's consistency conditions |22|j23[| . These are 
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T„o p ^„o p (k — q — 0) — - 



. qyiv(O) 
F2 



+7r o p (fc = o, 9 2 = m 2 ) = o , 

^o„(k 2 = m\,q = 0) = , 



(2.17) 

(2.18) 
(2.19) 



where ova/^O) is the well known ttN u-terni |24[] and F n the pion radiative decay constant. These three conditions 
yield only two relations amongst the constants, which are 



goo 



f]0c 



= ( — - ) {gia + g2a) 



(2.20) 
(2.21) 



Furthermore, we recall that the on mass shell 7r°p S-wave elastic scattering amplitude is parameterized according to 

PIS 



F = a+ + b+k • k , 



(2.22) 



with a + being the isospin even 7r°p scattering length. To obtain a third relation we may require that the amplitude 
Eqn. 17 reproduces F at k = 0, i.e., 



*n°p(k 



m^k 



0;q 2 



0)-«*r(l + £) 



This yields, 



goa 



2(^fg^ = 



An 1 



M 



(m c 



frrNN [jf 



AM 2 



AM 



J2 



(2.23) 



(2.24) 



We now may use <t w jv(0), F v , a + and g a NN as input to evaluate the vertex constants. The value of the quantity 
ovjv(O) is related to matrix elements of the operator m q qq in the proton, where m a stands for the mass of the proton 
quark constituents. This quantity can be calculated from the baryon spectrum p7| . To leading orders in the quark 

3/2 3/2 O 

masses one finds to order 0(m q ) that avjv(O) = 26 MeV. Including 0(m q ) contributions and estimate 0(m 2 ) 
ones yields, 



0Vjv(O) 



35 ±5 



1 



MeV 



y 



(2.25) 



with y being a measure of the strange quark content in the proton. An even higher a value is obtained from the 
isospin even 7rN scattering amplitude at the Cheng-Dashen point, which can be determined exp erime ntally using the 
low-energy theorem of current algebra. This gives a value ov jv (0) = 45 ± 8 MeV. Equation 2.25 suggests a value 
Pttn(0) = 35 MeV for y = 0, in keeping with the OZI rule pS] . With this value and taking a radiative decay constant 
= 93.5 MeV |2i|; isospin even 7r°p scatte ring l ength a^~= = (—0.010 ± O.OOS)™^ 1 p5| , p6| , and the erNN coupling 
g^^^/A-K = 8.03 |19| we can now solve Eqns. 2.2C , 2.21 , 2.24 to extract the values, 



g 0(T = 0.22 ± 0.04, g la = -0.82 ± 0.42, g 2a = -2.62 ± 0.7 



(2.26) 



With these constants Eqns. |2.16| predicts 6+ = (-0.068 ± 0.015)m~ 3 in close agree ment with the experimental value 
b+ xp — —(0.088 ± 0.014)m~ 3 quoted by Koch |25[| . Thus the amplitude Eqn. 2.16 on the mass shell reproduces the 
correct k • k dependence of the 7r°p elastic scatteri ng am plitude, Eqn. 2.22, also. 

It is difficult to ascer tain that the amplitude Eqn. 2.16| does reproduce the off mass shell behavior correctly. However 
we m ay co nfront Eqn. [2. 16 with predictions from other approaches. For example, taking the residue of the amplitude, 
Eqn. 2.16, at the a pole (q — k) 2 = m 2 ., one obtains an effective aim coupling 



7 e// 



V^Aq = mi, k 2 = mi, (q - k) 2 = m%) = (2.5 ± 0.9)m o 
a value to be compared with the well known estimate from soft pion physics |p2| 

™2 



V r (TW 7r(m w ,m 7r ,m .) = m a 



2F W 



1 



(2.8)m . . 



(2.27) 



(2.28) 
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Two-pion loops contribute to the tttt elastic scattering amplitude and appear as corrections to the leading contact 
term [122 30 1. As a further check we may evaluate these corrections in terms of an effective a- meson exchange in s, u 
and t-channels. By doing so it is straightforward to show that these corrections are 



rploop 



7T7T — >7T7r 



-V^(mi,mi,s)S ab 8 ca + 



m~ — u 



-V^{m^mi,u)5V + 



(2.29) 



where 6 l ' J is the Kroneker 5 and with i,j being pion isospin indices. The factor of 2 is due to the symmetry of the 
two ainr vertices. At threshold this expression amounts to about 20% of the contact term, in good agreement with 
the value of 25% reported by Gasser |}(| from improved low energy theorems. 

Now that the airir vertex parameters are well defined we may examine offshellness effects. To this aim we evaluate 
T w o p ^o p at the production threshold of p p — > p pir and compare with its on mass shell value. When both pion legs 



are on the mass shell the last term in Eqn. 2.16 reduces to 



T 



fjp^opik? = ml; q 2 = mi) 



■ 9aNN 



90a 



~i0.ll 



9<jNN 



At threshold of the pp — > ppir reaction, the momentum squared of the a meson is (q — k) 2 = (p' — p) 2 
the off mass shell t pole term in Eqn. 2.16 becomes, 



(2.30) 
-Mm-* and 



T 



(t) 



-igaNN 



goo 



-il.O 



m 2 + Mm v 

2 



/-(-Mm-) 



m„ 
g<jNN 
m„ 



gio 



Mm.* . 

2~ K9\a + 32<t) 



(2.31) 



which is a factor sa 9 larger compared to the on mass shell value Eqn. 2.3C . In Fig. 5 the amplitude T^Op^^Op is drawn 



as a solid line vs. q . The contributions from the s and u channel nucleon pole and a- meson pole t erms are drawn 
as dashed and dot-dashed curves, respectively. Off mass shell behavior of the s and u terms in Eqn. 2.1£ is given by 
the pion form factor, f^iq), the nucleon propagators and the momentum dependence of the nNN vertex. That of the 
t-channel is affected by the a meson form factor f a (q— k) and the Va-nn vertex function, Eqn. 2.15. Note that on mass 
shell both of these contributions are small, opposite in signs and cancel to large extent. As q z and (q — k) 2 become 
more negative both terms become negative and therefore add constructively to the conversion amplitude, giving rise 
to strongly enhanced off mass shell amplitude. We expect then that 7r exchange plays an important role in the 7r° 
production process. 



C. The T vp _> n0p amplitude 

The production of a neutral pion can occur also on an internal <5-meson line. To evaluate the amplitude for the 
conversion process r\p — > 7r°p we apply a similar procedure as above. The r\ meson couples to the nucleon isobar 
A* (1535 MeV) strongly |3lj and the main contribution to the conversion amplitude for r/p — > 7r°p is due to s and u 
isobar pole terms (see Fig. 6), hereafter we refer to as the resonance contribution. The other terms (graphs 6c-6e) 
furnish a background term. A t-channel would involve a vertex with 6, rj and ir legs defined to be, 



V5 V Tt{k,q) = m s 



9os 



Tgu 



5-32,5 



{k-qf 



-gs.8 



(2.32) 



Here, to account for the fact that the three legs are different the vertex is described in terms of four constants . Taking 
the sum of all the graphs in Fig. 6, we may write, 



6 



T 



rjp — >tt u p 

-ig v NN'girNN»f v {q)f7r(k)u(p') 

fr,(<l)U( k ) u (P 



1 



AfR-- v /s + *r/2 Mr- y/n + iT/2 

1 1 \ 1 

M 



u(p) 



AP 



u 



M 2 



u{p) 



igsNN , ^."L k 2 V5t,ir{k,q)u(p')u(p) . 
(q — ky — rug 



(2.33) 



Here Mr = 1535 MeV and T = 175 MeV stand for the mass and width of the isobar resonance, and the coupling 
constants are g V NN* — ^-^,g-KNN* = 0.8 [[H). 

We now apply the Adler's consistency conditions to the background term in Eqn. 33 to write the following relations 
amongst the unknown Srjir vertex constants (see Eqns. 17-19), 



505 



m s ovat-^a^O) 



F F 

± IT- 1 - Tj 



gsNNfs(o) 

9os = - [ — ) (.92,5 + g3s) 
\ms J 

2 



.905 



{gu + g3s) 



(2.34) 
(2.35) 
(2.36) 



where F„ is the r\ radiative constant; in the limit of an exact SU(3) symmetry F„ — F^. Likewise, the quantity 
v-nN^N is a cr-term, a quantity related to matrix elements of various quark mass qq in the nucleon, 



®ttN — >r)N (^) = —(p'\lLU + dd\p) 



(2.37) 



where to = (m u + md)/2 is the average of the u and d quark masses and t = (p' — p) 2 the transferred momentum 
squared. This term can be deduced from the kaon-nucleon (KN) cr-terms and the strange quark mass through 24 



OVA^ t; Ar(0) = 



m + m s 



(2.38) 



Taking the quark masses as m u = (5 ± 2) MeV, m d = (9 ± 3) MeV, m s = (175 ± 55) MeV, and the KN cr-tcrm 
cr^]v(0) ~ (200 ± 50) MeV , cr-term a^ ] N (0) ~ (140 ± 40) MeV , one obtains, 



(0) ~ (2± 1.3) MeV . 



(2.39) 



Now the partial decay width of the 5 meson into a ttt] pair is |2£|] F nri = 57 ± 11 MeV . By calculating this width 
using the expression Eqn. 32, one obtains a fourth relation amongst the Srjn vertex constants, 



gos 



m 2 



gu 



— .925 + .935 = 
v TO 5/ 



/r87r 

"W 



Resolving this last relation and the Adler's conditions, Eqns. 2.34 - 2.36] , gives 

g os = 0.05 ± 0.035 , g is = -5.95 ± 1.5 , g 2 s = -3.46 ± 0.1 , g 3S = 3.30 ± 0.2 



(2.40) 



(2.41) 



With these constants one finds that at threshold for the pp — > ppir reaction, the conversion amplitude is T r]p ^^a v w 
0.98 /to. The t-pole contribution amounts to only T^J^^o ~ 0.04 fm We may thus conclude that a S meson pole in 
a t-channel contributes very little to T vp ^^o p and should play a minor role in the ir° production process. 



D. The T pp _^ 0p amplitude 

We follow van Kolck et al. M| and limit the pun vertex to the form, 
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(2.42) 



where the coupling constant g pul7r ~ —10, a value fairly well established. With this expression, the contribution from 
a puj exchange mechanism to the Tp^^^o pj is 



1 



(2.43) 



1 



(q x p + ipq ■ a — iaq ■ p) 



1 



_E + M ^ r ' E' + M 

The s and u nuclcon pole terms can be calculated from Eqns. 10 and 14. 



(q x p' + ip'q ■ a — iaq ■ p') 



III. S-WAVE AMPLITUDES AND CROSS SECTION 



In this section we write the primary production amplitude in a form suitable for numerical integration. We call 

Pj 



n, 



Ej + M 



(3.1) 



where pj and Ej are the three momentum and total energy of the j-th nucleon. For the incoming particles in the 
center of mass system (CM) III = II2 = II. The energy available in the CM system is 



Q = \fs — 2M — m n 



(3.2) 



where s = (pi +P2) 2 is the total energy squared. We shall calculate the amplitudes and cross sections as functions of 
the variable Q. We also define 



-M (m* + Q) 



a = 1 
6=1- 
R 



Q 



Q 



1 



M 



1 + — + 

2M V 2M 4M 2 J ' 



1 



1 



M R - M - - Q + iT/2 Mr - M + m n + Q + iT/2 
Following the discussion in the previous section we write the primary amplitude as 

M^ = M-x + Mr, + M a + M p + M u + M s , 



(3.3) 
(3.4) 

(3.5) 



where Mb stands for a partial production amplitude representing the contribution from the exchange of a boson B. 
We substitute equivalent two-component free spinor matrix elements (Table A6.1 of Ref. 26|) in Eqns. 2.11-2.13 for 
s-channels and the appropriate ones for u-channels, to write the scalar and vector meson exchange amplitudes as, 



Ms 



J-kNN 9s nn fit \ m 7r + 2Q M 

,2 _ „2 J SKI) 



and 



My = i 



TUn m s — q*- 



f-xNN 9vNN 



M + Q 2M + m v 



II • (01 - a 2 ) , 



m n my — q 2 



f v (q) [X II • (01 - a 2 ) + 1Y II • a x x <r 2 ] 



(3.6) 



(3.7) 



where we have used the notation 

X = 



2ll+ 3 -^]+ Kv (l+ m ^ + Q 



2M 



3™^ + Q nhr + Q . .. 

f K V — r^r — I 1 



2M 



1 - 



Y 



M 

-2-2(1 + k v ) 



2M V 2 ( 2M + ™>„) 
( StUtt + Q \ ( + Q 



2{2M + m n ) 



(i-K V n-n) 



V M 



2M 



rn n 



2(2M + m* 



(3.8) 
(3.9) 
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The puj exchange mechanism contribute to the production amplitude a term, 

M jt) = 2 ^L_ /2(g)2 _^L f M)gp ^ {p . U )U u\ X (J2 • (3.10) 

' nip — q z H m^j — q z m w 

To avoid double counting we add this to M p (but not to M u ) only. For both of the it and 77 exchange amplitudes, 
there is a t-channel pole in addition to s and u-channel nucleon pole terms. The sum of all three gives, 

/2M\ ( 1 
\m n J \ml - q 2 



Mir = -if„NN —5 o 

\m„ J \mi - q l ) 



Mr] = -ifxNN ( ~Z7T — 2 ) Mi) 



(3.11) 



\ m n / \ m 2 — q- 

9riNN*9irNN*R + firNNfnNN ( (-) + gSNN ( ^2^2 I ^V^(^ 2 = m Tr' I 2 ) 



(3.12) 



Finally, the total cross section is calculated from the expression 

M 4 



o = / ^ fP4 ^P, 1 ZM (») | 2(J4( } (3 13) 

16(27r) B v/(«)p x J E 3 £4 E v 1 1 ^ Pih V ; 

where Z is the three-body FSI correction factor of Refs. (20) to be specified below. 

IV. PREDICTIONS AND COMPARISON WITH DATA 

We apply now the model presented in the previous sections to calculate the total cross section for the pp — * ppw 
reaction at energies close to threshold. We first consider the relative importance of the various exchange contributions. 
To this aim we draw in Fig. 7 the primary production amplitude M^ m ) along with the partial exchange amplitudes 
of Eqns. 46-51, vs. the energy available in the center of mass (CM) system. The main contribution is due to 
pion exchange with ratios : M p : : M a : : M$ 138 : 42 : 8 : 6 : 5 : 0.6. Next important to the 

pion is contribution from the p meson {M { p s+U) « 30% and M { p t] » 4.5% of M^). In comparison with M n other 
contributions are significantly smaller. Nonetheless, they influence the cross section strongly through interference. 
They all have a common relative phase and add constructively. The p, having an opposite sign counteracts to balance 
their effects. 

The ratios quoted above differ considerably from those predicted by Horowitz et al. ||. Particularly, the s and 
u-channel nucleon pole terms for scalar meson exchanges have opposite signs, thus suppressing the a exchange con- 
tribution and practically eliminating that from the 6 meson. We note however, that as in Ref. j^] the nucleon pole 
term contributions from n and a exchanges add constructively. 

We want to emphasize at this stage that, our transition operator M^ m ) for the production process, involves contri- 
butions from connected diagrams (Fig. 3) only, and accounts for all relativistic and crossing symmetry effects. Taking 
u-channel contributions only, one obtains far more important effects than with both of the s and u-channels. The 
lesson to be learned here is that many small terms added up coherently can explain a seemingly large discrepancy, at 
a time when other small terms which counteract to balance their effects are disregarded. Under these circumstances it 
seems essential to preserve crossing symmetry at all stages of the calculations and treat all contributions on an equal 
basis. 

It is rather surprising that the a contribution is not as important as predicted in Ref. ||. Yet, it is to be indicated 
that in Ref. ||, only contributions from negative-energy intermediate states are included explicitly to the pion pro- 
duction operator via a Z-graph, while contributions from positive-energy nucleon intermediate state are presumably 
contained in the distorted waves describing the initial and final two-nucleon states. Thus, different approximations 
are used in the calculations of the direct and Z-graph from either s or u- channels and it is not clear how a delicate 
balance between these contributions is maintained throughout the calculations. 

In Fig. 8 we draw predictions for the total pp — ► ppn cross section as obtained for three different value of the ttN 
cr-term. The small dashed curve shows results obtained without the tr-meson pole term included. Clearly, the overall 
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contribution from nucleon pole terms only does not provide the enhancement required to resolve discrepancy between 
the calculated cross section and data. Albeit, the contribution from an isoscalar-scalar a- meson pole term dominates 
the production cross section and as shall demonstrate below, with a ovat(O) = 35 MeV, the a pole term provides the 
enhancement required to explain cross section data near threshold. 

The strong dependence of our predictions on the ttN er-term deserve a comment. Though still not very well known, 
it seems more likely that the value of this quantity falls in the range ovat(O) = 35 — 45 MeV ^7j. In the calculations 
presented in the present work we have not included contributions from A(1232 MeV) isobar excitations. It remains 
still to be verified whether the effects of such term would bring the calculated cross section to agree with data even 
with a larger value of ovat(O). 

V. FSI AND COMPARISON WITH DATA 

In comparison with data, the predicted cross section in Fig. 8 varies very fast with energy due to phase space factor. 
The primary production amplitude (see Fig. 7) is practically constant near threshold and therefore contributes very 
little to the energy dependence. To cure this deficiency of our predictions and allow for comparison with data to be 
made we must account for ISI and FSI effects. In what follows we treat ISI and FSI using two different approxima- 
tions, which as we shall demonstrate below both yield similar cross sections. 

Approximation! 

In this approximation the transition operator M^ m > is treated as an effective operator acting on nucleon wave func- 
tions. These are calculated using a phenomenological NN potential. In this approximation we neglect interactions of 
the 7T° produced with the outgoing nucleons. This is a standard and usual procedure in the literature [fUTsll, reducing 
a three-body problem into effectively a two-body process. To be consistent with the OBE picture applied in the 
present work, we use the OBE NN potential of Machleidt jl!|; potential parameter set C of his table A. 2. The initial 
and final two-nucleon radial wave functions are calculated in momentum space from half off shell R matrix, following 
exactly the procedure of Ref . |l9[ . The evaluation of the transition amplitude and cross section is performed as in 
Rcf. We have verified that the calculated wave functions reproduce well the known experimental phase shifts 

& 

Approximation! 1 

Here the S-wave production amplitude is assumed to factorize into a primary production amplitude M^ ln \ a P-wave 
ISI factor and and an S-wave FSI factor [^o) , 

T 2 ^ 3 - (^P\Mti\<;.! (2) ) « Z 33 Mtlx 2 2 , (5.1) 

where ^ el \{2) and ^^(3) represent the initial two-body and final three-body wave functions. The ISI correction 
factor is taken to be X22 ~ |1 + sin(53 Po ) exp i(53 Po )\. Here <$3p is the pp P-wave phase shift. At threshold of the 
pp — > ppir reaction, 8zp w —10 so tha ISI corrections to the cross section are bound to w ±30%. For a three 
body process as in our case, the FSI correction factor is identified pp|] with the elastic scattering amplitude (on mass 
shell) for the process nNN — > irN N (three particles in to three particles out) and has the structure of the Faddeev 
decomposition of the t- matrix for 3 — > 3 transition. An important property of this approximation is that the different 
two body interactions among the out going particles contribute coherently. Although the meson-nucleon interactions 
are weak with respect to the NN interaction, they can still be influential through interference. This approximation 
was discussed in length elsewhere [^0| and we skip further details here. In the analysis presented below the FSI factor 
is estimated from 7rN and NN elastic S-wave scattering phase shifts |34|,[35|]. The S-wave NN phase shift is obtained 
from the effective range expansion which includes Coulomb interaction between the two protons. We have used the 
scattering length a pp = —7.82 fm and an effective range r pp = 2.7 fm of Ref_j3J]. The Sll and S13 7rN scattering 
lengths are taken to be a\ = 0.173 m^ 1 and 03 = —0.101 m~ x , respectively |35| " 

To compare the results from the two approximations we draw in Fig. 9 the partial cross sections corresponding 
to production via a tr meson t pole mechanism. For neutral pion production, the cross sections calculated using 
the two approximations are practically identical. We note though that Approximation I does not account for final 
meson-nucleon interactions. To see how influential the 7rN FSI can be, we draw in Fig. 10 the cross section corrected 
for FSI interactions assuming pure 1=^ and I=| interactions for the outgoing 7r°N pairs. The energy dependence 
of the cross section differ significantly for the two channels. But, when the 7rN interactions are taken in the proper 
isospin combination, their overall contribution to the FSI factor almost cancels out, leading to an energy dependence 
practically identical with the one obtained with charged pp FSI (small dash curve) only. We want to stress here 
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that this is an accidental consequence of the fact that the 7r°p interaction over the allowed 1=1/2 and 1=3/2 isospin 
channels averages to zero, making contributions from diagrams with 7T°p interactions ineffective. This may not be the 
case for charged pions and r] meson production ppfl . 

VI. SUMMARY AND CONCLUSIONS 

We have calculated S-wave pion production in pp — ► ppir using a covariant OBE model, where a boson B created 
on one of the incoming protons is converted into a neutral pion on the second. The amplitude for the conversion 
process BN — > ir°N are taken to be the sum of s and u pole terms and when allowed a meson pole term in a t-channel. 
To be consistent with the OBE picture of the NN interaction we have considered contributions from all of the n, r), 
a, p, to and 6 mesons. Both the scale and energy dependence of the cross section are reproduced rather well. 

Based on a covariant quantum field theory, the transition operator must be fully covariant and as such it includes 
contributions from irreducible connected diagrams only, each involving large momentum transfer. The operator M^ m \ 
Eqn. 48, and likewise the calculated amplitude preserve crossing symmetry and relativistic covariance throughout the 
calculations. Furthermore, terms which involve three meson vertices are written in a general form. As an example, 
the Lagrangian corresponding to a arnr vertex, in addition to non derivative terms, includes derivative a(dirdn) and 
(<9er<97r)7r interaction terms. 

The off-shell properties of the various partial production amplitudes are taken into account explicitly. For vector 
and scalar mesons, offshellness is contained in the nucleon form factor and has a marginal influence. A pseudovector 
coupling for the pseudoscalar particles, as assumed in this work, introduces an extra factor of q to the amplitude. 
Consequently for the tt and r), each of the s and u nucleon pole terms at q 2 — — 3.3/to -2 becomes a factor of « 7 
higher in comparison with its value on mass shell. To account for the off mass shell behavior of the cr-meson pole 
term, we have used Adler's consistency conditions and the isospin even 7r°p scattering length to write the cnrir vertex 
in a rather general form. Taking all of the s, u and t-channel contributions into account, we have found that the 
ir°p — » ir°p conversion amplitude off the mass shell is more than a factor « 15 higher in comparison with its on mass 
shell value. Consequently, the tt exchange contribution to the amplitude for the pp — > ppir reaction dominates the 
production process and brings the calculated cross section to agree with data. 

We have found that the contribution from a cr-meson pole in a t-channel dominates pion production at threshold. 
A term as such accounts effectively for two-pion exchange contributions. The <77T7r vertex used to calculate this 
contribution depends rather strongly on the 7rN cr-term. However, even if we use a value, 0v./v(O) = 25 MeV, which 
is unrealistically too low, a t-channel a pole term amounts to more than a factor of two higher a contribution than 
that of the s and u-channel nucleon poles. It seems quite impossible to disregard such an important contribution 
to the production process. Neither the overall contribution from all meson exchanges, nor any of the various meson 
exchange contributions, from s and u nucleon pole terms only, can reproduce data for this process. In fact, various 
contributions of a similar size are found to interfere destructively and including all of these in a consistent manner, 
yields a production amplitude far below what would be required to explain data. This conclusions seems unavoidable 
and independent on the approximation used to account for ISI and FSI corrections. 

A meson production in NN collisions necessarily involves large momentum transfer. Thus, two-pion exchanges 
which describe short range interactions are expected to play an important role. In a traditional covariant OBE model 
such contributions are represented by an effective scalar cr-meson pole term. In xPT calculations these should result 
from expanding the existing calculations to include one loop two-pion diagrams. Recently, we have extended the 
existing xPT calculations |Tj|-|l3| taking into account all tree and loop diagrams up to chiral order D=2. With these 
one loop contributions added, it is found that in analogy with a dominant a meson pole term, there are substantial 
contributions from isoscalar-scalar two-pion t-channel exchanges . 

In summary, two-pion exchanges play an important role in the production of ir° in pp —> ir°pp. By taking these 
into account, we have demonstrated that a covariant OBE model can explain the existing near threshold neutral pion 
production data . 
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FIG. 1. The primary production mechanism for the NN — > NN-k reaction. In this figure and following figures, solid lines 
represent nucleons and broken lines mesons. 
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FIG. 2. Graphs contributing to the amplitude of the conversion process BN — ► ir°N. Graphs a and b stand for nucleon pole 
in s and u-channel; the graph c represents a B' meson pole in a t-channel. 
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FIG. 3. The production amplitude. The mechanisms of graphs a and b represent pion production on external nucleon lines. 
Graph c represents production on an internal meson line. From parity and isospin conservations the latter mechanism is limited 
to production from a and 5 meson lines only. 
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FIG. 4. Terms contributing to the it p — ► 7r p conversion amplitude. Graph c describes an effective isoscalar a-meson pole 
in a t-channel. 
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FIG. 6. Terms contributing to the conversion amplitude for the rjp — ► 7r p. Because of the strong coupling between the 
^-meson and the iV* (1535 MeV) nucleon isobar one should allows for contributions from nucleon and nucleon isobar intermediate 
states. Graph e describes 5-meson pole in a t-channel. 
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FIG. 7. The primary production amplitude for the pp — > ppir reaction vs. Q, the energy available in the overall center of 
mass system. Different exchange contributions to M' <n ' are drawn separately. 
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FIG. 8. Predictions for the total cross section of the pp ppir° reaction for different values of cr^jv(O). Predictions with s 
and u nucleon pole terms only are shown by the dotted curve. All curves are not corrected for FSI. 
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FIG. 9. ISI and FSI corrections. Energy integrated cross sections assuming production via a a meson t pole mechanism only 
with Approximation I (solid line) and Approximation II (dot-dashed line) 
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FIG. 10. FSI corrections. Integrated energy cross sections calculated with the assumption that the interacting 7rN pair is 
scattered in isospin 1=^ (large dashed curve) and I =| (dot-dashed curve). The solid line is that obtained with the ir N 
interactions taken in the appropriate isospin combinations. Predictions which account for the pp FSI only (small dashed curve) 
are nearly identical with the solid curve. All predictions are calculated with av^O) = 35 MeV. The data points are taken from 
Refs. H 
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